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ABSTRACT 

The instability of large amplitude Gortler vortices in a growing boundary 
layer is discussed in the fully nonlinear regime. It is shown that a three- 
dimensional breakdown to a flow with wavy vortex boundaries similar to that 
which occurs in the Taylor vortex problem takes place. However, the instabil- 
ity is confined to the thin shear layers which were shown by Hall and Lakin 
(1987) to trap the region of vortex activity. The disturbance eigenfunctions 
decay exponentially away from the center of these layers so that the upper and 
lower shear layers can support independent modes of instability. The struc- 
ture of the instability, in particular its location and speed of downstream 
propagation, is found to be entirely consistent with recent experimental re- 
sults. Furthermore, it is shown that the upper and lower layers support wavy 
vortex instabilities with quite different frequencies. This result is again 
consistent with the available experimental observations. 


This research was supported under the National Aeronautics and Space Admini- 
stration under NASA Contract No. NAS1-18107 while the authors were in resi- 
dence at the Institute for Computer Applications in Science and Engineering 
(ICASE), NASA Langley Research Center, Hampton, VA 23665. 
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1. INTRODUCTION 


Our concern is with the nature of the three-dimensional breakdown of 
steady, spanwise periodic large amplitude Gortler vortices. It is known from 
the experiments of Bippes and Gortler (1972) and Aihara and Koyama (1981) that 
this breakdown leads to a time-periodic flow with wavy vortex boundaries simi- 
lar to those which occur in the Taylor problem. More recently Koharaa (1987) 
has investigated vortex instabilities in boundary layer flows over a laminar 
flow wing and found a secondary instability of Gortler vortices localized at 
the top of the region of vortex activity. Furthermore Kohama found that the 
instability propagated downstream with a speed which approached the free 
stream speed as it developed in the downstream direction. The onset of this 
time-dependent motion was found in all of the above experiments to be ulti- 
mately followed by transition to turbulence. At this stage more obvious dif- 
ferences between the Gortler and Taylor problems emerge so that, for example, 
the rich bifurcation structure of the Taylor problem is apparently not carried 
over to the Gortler problem* 

In fact even in the linear regime the apparent similarities between 

Taylor and Gortler vortices are perhaps misleading since it is known from the 

work of Hall (1982a, b, 1983, 1984) that nonparallel effects in the Gortler 

problem cannot in general be ignored. Indeed it was shown by Hall (1983) that 

the inconsistences of the various parallel-flow theories of for example 
»• «• 

Gortler (1940), Hammerlin (1956) and later authors are a direct consequence of 
the parallel-flow approximation. The only regime where this difficulty with 
the parallel-flow theories does not occur is at small vortex wavelengths where 
the effect of boundary layer growth on the vortices becomes less important. 
However, the most surprising feature of the nonparallel theory of Hall (1983) 
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is that the concept of a unique neutral wave is not tenable in the Gortler 
problem since the downstream position where a vortex begins to grow is a func- 
tion of the location and initial form of the imposed disturbance. A more sig- 
nificant consequence of the nonparallel theory is that the concept of a unique 
growth rate at a given downstream location is also not tenable, this result 
makes transition prediction by empirical methods such as the e n rule not 
possible. Nevertheless, there is much work still being done in the context of 
"parallel-flow" Gortler vortices. 

The nonuniqueness properties of solutions of the correct zeroth order 
approximations to the linear Gortler vortex equations were shown by Hall 
(1986) to occur in the corresponding nonlinear problem at 0(1) wavenum- 
bers. In the latter paper, the development of finite amplitude Gortler vor- 
tices was investigated numerically using a finite difference discretization in 
the normal and chordwise directions together with a Fourier expansion in the 
spanwise direction. It was found that as the vortices move downstream the 
disturbance energy of the flow becomes concentrated in the fundamental and the 
mean flow correction. This is entirely consistent with the weakly nonlinear 
theory of Hall (1982b) which is appropriate to small wavelength vortices. 
However in a growing boundary layer, the "local 11 wavenumber of a fixed wave- 
length disturbance grows like the displacement thickness of the boundary layer 
so that in most flows any vortex will eventually enter the small wavelength 
regime where locally the asymptotic analysis of Hall (1982a, b) apply. 

The surprising feature found by Hall (1987) in a numerical simulation of 
»« 

nonlinear Gortler vortices is that even at relatively low wavenumbers the dis- 
turbance organizes itself so that almost all of the energy is in the funda- 
mental and mean flow correction. This is precisely the situation found by 



-3- 


Hall (1982b) in an asymptotic investigation of nonlinear Gortler vortices 
where it was shown that the Stuart-Watson type of description of the onset of 
finite amplitude motion is replaced by a "mean-f ield" interaction. 

As is the case with all weakly nonlinear stability calculations the work 
of Hall (1982b) is restricted to a neighborhood of the position where a given 
disturbance is neutrally stable. However, it can be inferred from that calcu- 
lation that a vortex of nondimensional wavenumber e * where 0 < e « 1 
reinforces the basic flow at zeroth order at a distance 0(e) downstream of 
the neutral point. This result was recently developed by Hall and Lakin 
(1987), hereafter referred to as HL, to give an asymptotic description of 
fully nonlinear Gortler vortices at 0(1) distances beyond the neutral 
point. It is the instability of this type of vortex which we will investigate 
in this paper. However, before discussing the nature of this instability we 

need to point out the salient properties of the HL calculations. 

-1 

Consider then a Gortler vortex of wavenumber e developing in a 

•• -.4 

boundary layer flow with Gortler number of 0(e ). This choice of small 

wavelength vortices is not as restrictive as it might first appear since, as 

explained above, this regime is always approached by a fixed wavelength vortex 

in a growing boundary layer. Suppose further that the flow is neutrally 

stable at the downstream location x = x , then HL showed that for x > x 

n’ n 

the flowfield splits up as shown in Figure 1. The vortex activity is confined 
to region I and decays exponentially to zero in the thin shear layers IIa,b. 
In regions IIIa,b there is no vortex activity and the mean flow satisfies the 
boundary layer equations. However, in region I the mean flow is determined as 
a solvability condition on the equations for the fundamental. In fact, the 
mean flow adjusts itself so as to make the fundamental and all the higher har- 
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The mean flow equations then determine the vortex 
velocity field in I so that there is a complete reversal of the usual roles of 
the mean flow and harmonic equations compared to say the situation in flows 
where nonlinearity can be described by the Stuart-Watson method. The shear 
layers located at y^ and y2 change position as they move downstream; their 
positions are determined from the solution of a double free boundary problem 
associated with the boundary layer equations. However, in flows where the 
local Gortler number increases faster than the fourth power of the local wave- 
number HL showed that y^ migrates to the wall whilst y2 moves to or beyond 
the edge of the boundary layer. The mean downstream velocity components in 
the layers IIa,b then approach the free stream speed and zero respectively; 
this has fundamental implications for the time dependent structure of the 
breakdown of this flow. 

We shall seek secondary instabilities of the flow in layers IIa,b; more 
precisely we superimpose spanwise periodic travelling waves on the flow in 
these layers and see how they develop. These perturbations are y radians 
out of phase with the fundamental so the secondary instability if it occurs 
will produce locally wavy vortex boundaries in IIa,b. It is of course not 
obvious that, should wavy vortices occur, the regions IIa,b should be 
particularly susceptible to these modes. In order to see why this is the 
case, we consider the model equation 

ay 

together with the condition + 0 , |y| °°. 


I 



-5- 


In fact, this equation essentially governs the nonlinear growth of time- 
dependent Gortler vortices in curved channel flows. Here X is a 
parameter, t denotes time, x is the distance around the channel, and y is 
the distance from the center of the internal viscous layer where the vortices 
initially develop. A finite amplitude solution of (1.1) representing a 
steady, x independent vortex i|> * satisfies 



" 4“ + X},| v 



0 . 2 ) 


and the instability of this flow to a travelling vortex like perturbation can 
be calculated by setting 

* = * v + r(x,y,t) 


where is a real function of x, y and t. If we now look at dis- 

turbances with ty' = Real(e* kx+at <)>(y)) then we find that 4>(y) satisfies 

a 2 2 2 

( — j “ J~ + X - iky - o)4> = 3i}> v <j> , <t> -*■ 0, |y| > ». (1.3) 

3y 

which determines an eigenrelation a = a(k). However if the wave is y 
radians out of phase in the spanwise direction with the fundamental vortex 
then (1.3) becomes 

a 2 v 2 2 

( — j - -t— + \ - iky - a )<j> = iji <|> , <(>-»- 0 , | y | 00 . (1.4) 

9y v 

The latter eigenvalue problem was studied by Shaw (1985) who found that 
unstable modes occur for A = 0(1). (In contrast to this result (1.3) leads 
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to stable disturbances.) Here X plays the role of Gortler number so 
increasing X is equivalent to increasing x for the boundary layer prob- 
lem. Of more relevance to our problem is the solution of (1.4) for X » 1. 

In this case \|» v develops a triple layer structure with a core region 

2 1/2 -1 / f > 

with i|> » (X - y /4) ' trapped between shear layers of thickness X 

1/2 

at y = ±(4X) . These shear layers correspond to IIa,b in the HL calcula- 

tion. An examination of (1.4) in this limit shows that any eigenvalues must 
now concentrate in IIa,b because in the core <j> now satisfies 

3 2 

( — y — iky — <7 )<p — 0. 

9y 

This is an Airy equation so that <(>•*■<» when y + 2/T or y -> -2/X 
and matching with the corresponding solution in the shear layers cannot be 
achieved. Thus the eigenfunctions must now concentrate in the shear layers 
and decay away from the centers of the layers. 

The above structure for the model equation is sufficiently close to the 
boundary layer problem for it to be applicable there. Thus, after formulating 
our instability equations in Section 2 we will in Section 3 investigate the 
instability of IIa,b to wavy vortex modes. In Section 4 we present the re- 
sults of our calculation whilst in Section 5 we draw some conclusions. 

2. FORMULATION OF THE PROBLEM 

The flow under consideration is that described by HL. We consider the 
flow of an incompressible, viscous fluid of kinematic viscosity v and 


i 


density p, 


over a wall of variable concave curvature 


a *x(X/L). Here 
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a Is a typical radius of curvature, X denotes distance along the wall and 
L is a typical length scale along the wall. The Reynolds number for the 

flow, Rg, is defined by 



( 2 . 1 ) 


where Uq is a typical flow velocity. A curvature parameter, 6 c> 
by 

« = - . 

c a 


is defined 
( 2 . 2 ) 


We are interested in the limit 
by 


Rjy -*■ <*> with the Gortler number G, defined 


G = 2r£ /2 <S ( 


( 2 . 3 ) 


held fixed. We denote time by T and (X,Y,Z) are taken to be the co- 
ordinates along the wall, normal to the wall and in the spanwise direction re- 
spectively. If (U,V,W) denotes the corresponding velocity vector we define 
dimensionless co-ordinates (x,y,z) and velocity (u,v,w) by 

(x,y,z) - L _1 (X,YR^ /2 ,ZrJ /2 ), 

and 

(u,v,w) = Uq 1 (U , VRg /2 ,WRg^ 2 ) . 


Our analysis is restricted to flows with u * 1 when y -► 00 and the 
pressure P is written in the form 
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The continuity equation and non-dimensional unsteady Navier-Stokes equa- 
tions for the flow take the form 


3u 3v 3w 

3x 3y 3z 


0, 


2 2 

3u , 3u 3u 3u 3 u 3 u 

3t dx 3y 3z ^ 2 ^ 2 


3y 3 z 


3v . 3v , 3v , 3 v 1 r 2 3p 3“v , 3~v 

_ +u __+v^ 7 +w^ r =-^G X u " ^ + 


3 v 
3x 


3 v 
8y 


3 v 
3 z 


3y 


2 » 2 

— 2 + — 1 ’ 
3y 3 z 


(2.4a,b,c ,d) 


2 2 

3w , 3w . 3w . 3w 3p 3 w 3 w 

FF +u 31T +v ^7 +w TI = -77 + ^T + ^7- 


Here the non-dimensional time variable t is given by t - UqL an< ^ terms 

of order Rl 1 ^ 2 have been neglected. HL obtained a steady solution of 
E 


(2.4) which satisfies 


u = v = w = 0, y = 0, 
u -»• 1, y <*>. 


(2. 5a, b) 


This was an asymptotic solution valid in the limit of small vortex wave- 
length. The Gortler number G, defined in (2.3), is expanded in the form 

G = GqS ^+Gj£ 2 +**», (2.6) 

where e~* is the non-dimensional wavenumber of the vortices. In the main 
part of the boundary layer u(x,y), the zeroth order, z-independent part of 
the downstream velocity component, satisfies 
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_ - a u 
G o xu 97 


1. 


which reflects the fact that here the mean flow is driven by finite amplitude 
vortices. The velocity field of the (smaller) vortices is then found by con- 
sidering the z-independent part of the equations of motion in I. This calcu- 
lation shows that the vortices are trapped between yj and y 2 and formally 
decay to zero exponentially in IIa,b. Above y 2 and below y^ there is no 
z-dependence to the flow and it is obtained by solving the boundary layer 
equations with jump conditions at y^, y 2 » The shear layers IIa,b correspond 
to the X layers for the model problem and we now examine the flow 
structure in Ila. In fact our analysis is equally applicable to region lib so 
our theory determines the stability of both shear layers which trap the 
vortices. In order to investigate the instability of the boundary layer in 
this region we consider perturbations to the steady, basic flow satisfying 
(2.4) and (2.5) in region Ila. 


3. THE ASYMPTOTIC STRUCTURE OF THE WAVY MODES 

2/3 

It was shown by HL that layers IIa,b are of thickness e so that in 

Ila we write 

(y-y 2 ) 

5 2JT~ ’ 

£ 


where y^Cx) is the location of the layer Ila. 

a _ y 2 9 9 u „ 1 8 

by 97 ~TTS w 97 by J 27 I H * 

in Ila, satisfying (2,4) and (2,5), are given by 


3 

Thus, in Ila we replace 
The basic, steady expansions 
equation (3.10) in HL and can 


be written in the form 
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u = Ug = Uq + e 2 ^ 3 Uj + ••• + e^ 3 cos(^)(Uqj + + •••) 

+ e 8/3 cos(— )(U + e 2/3 U +...) + ••• , 

e 02 12 

= , 2/3 = -2/3 ,Zw„ 2/3 . 

v = v g = v Q + e Vj + ••• + e cos(— )(V Q1 + e + •••) 

(3. la,b,c,d) 

+ e 2/3 cos(|^-)(V 02 + e 2/3 V 12 +•••)+••• , 
w = Wg= e 1/,3 sin(-^)(W 01 + e^Wj^ •••) + e 1 ^ 3 sin(-|^-) (W Q2 + e 2/ ^ 3 W 12 + ***^ + ***» 

P = P B = p 0 + e2 ^ P 1 + *** + e 4/3c ° S (|K P oi + ^ 2/3p n + •**) 

+ e 2/3 cos(|^)(P 02 + e 2/3 P 12 + •••) + •••, 


where the coefficients are functions of x and £• Note that the coeffi- 
cients u jk> V jk» Wjk> Pjiq in the above formulation correspond to 2Ujk> 

2iWjk, ^^jk respectively, j = 0,1,2,***, k = 1,2,3,**«, in the 

formulation of HL. The equations obtained from substitution of the basic ex- 
pansions (3.1) into the continuity and Navier-Stokes equations (2.4) and 

equating coefficients of like powers of e are given in detail by HL. The 
equation obtained which determines Vqj, from (3.1b), is 


where 


3^v 

_2i + g 5 v - 01 

35 2 *¥ v oi 


- s 2 fV 01 , 


i X'(a+2y„) 

«i« ’ = 


3/2 


- b)/3, 


3x^ G q x 


(3.2) 
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and 

/a+2y7 

g 2 ( x ) ■ 3 • (3.3a,b) 

Here f(x) is a function which can only be determined at higher order, a(x) 
and b(x) are arbitrary functions of x arising from the solution in region 
I (see (3.7) in HL), and a dash denotes a derivative with respect to x. 

The results of Davey, DiPrima and Stuart (1968) show that the Taylor- 
vortex flow is unstable against perturbations differing in phase from the 
fundamental component of the steady vortex flow by ^ . After instability, 
the new flow has wavy surfaces traveling in the azimuth separating neighbour- 
ing vortices. This suggests that, since we are seeking a secondary instabil- 
ity that will produce locally wavy vortex boundaries in IIa,b, we must con- 
sider a out-of-phase, time-dependent perturbation to the basic flow in 

Ila ,b. 

We have also considered the case of an in-phase perturbation and found 
the resulting problem more complicated than the present one. We anticipate 
that there are no unstable solutions for this case and just proceed with the 
present situation. 

Hence, we seek solutions with out-of-phase perturbations proportional to 

E = exp(-^y J X K(x)dx - ^) , (3.4) 

8 8 

where the wavenumber K expands as 

K = K Q + e 2/3 K 1 + ••• , 


(3.5) 
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and ft is the constant frequency. The length scale and time scale in (3.4) 

2 

9 ^ 3 9 9 

are chosen, from Hall (1982a), so that — s- ~ - 5 — and -r— + u -r— = 0 in 

OX OL OX 

O Z 

the shear layer, the latter scaling ensures that the waves travel downstream 
with the speed of the fluid in the shear layer. 

We find that the appropriate expansions in Ila take the form 

4/3 z 2/3 

u » Ug + ( 6 (e sin(— )E (uqj + e Ujj + •••) + •••) + ••• + c.c.}, 

v * Vg + ( 6 (e 2 ^sin(^-)E(vQj + e^^v^ + •••) + •••) + ••• + c.c.}, 

(3.6a,b,c,d) 

w = Wg+ { 6 (e *^cos(^-)E(Wqj+ e^^w^+ •••) + e ^^E(w ffl Q+ e ^^ w m i + •••) + •••) 

+ . . . + c.c .} , 

-4/3 z 2/3 

p = P B + { 6 (e ' sin(— )E( p q1 + e ' Pjj + •••) + •••) + ••• + c.c.}, 

where 6 is a small amplitude and c.c. denotes complex conjugate. Note 
that there are only mean flow (independent of z) correction terms occurring 
in the expansion for w as a result of the perturbation being y out of 
phase with the fundamental vortex. For the case of an in-phase perturbation 
mean flow correction terms also occur in the corresponding expansions for u, 
v, and p and the resulting eigenvalue problem is more complex. 

The coefficents are functions of x and £ and are determined by sub- 
stitution of the expansions (3.6) into the unsteady equations of motion 
(2.4). The zeroth order equations obtained from equating coefficients of 
6 sin(^-)E in (2.4a,b,c) and 6 cos(-|-)E in (2.4d) are 
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3 v, 


01 


- w ni " °» 


35 "01 


3 U 1 

•i(lu 0 i ♦ iK 0 = 0 u 01 ♦ v 01 Ti u 01> 


■“ v oi + 1K oVoi ’ ' VVoi ■ v or 


(3.7a,b 


0 " -Pm - w , 


01 01 * 


Consistency of (3.7b, c) requires that 


a = V 0 , 


which shows that the waves move downstream with mean speed Uq, i.e. 
speed of the mean part of the basic flow in Ila. Thus (3*7) reduces to 


9 v, 


01 


~ w m 


3 £ 01 


3Uj 

u oi + W V 01 ■ °- 


V 01 + fyVoi * °- 


”oi + p oi ' °- 


(3.9a,b 


= * u i 

(Note from HL that G xu n r; — = 1). From equating coefficients of 

0 U 

z 

6sin(~)E in (2.4b, c) at the next order we obtain the equations 


>c,d) 

(3.8) 
, the 

* c j d ) 
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_ 3u^ ^ u 2 

iK 0 U l U 01 + iK l U 0 U 01 + V 11 + W V 01 ” w m0 u 01 


3 2 u 


01 


35 


‘11 


8p 01 8 V 01 

1K 0 U 1 V 01 + iK l U 0 V 01 ■ W m0 V 01 = ' Wll " G 0 XU 1 U 01 - — + ^2" " v ir 

(3. 10a,b) 


It is useful here to note that from (3.12) of HL 


/ a+2y. 


> u, = 

•'GqX /a+2y 2 


/G 


0* 


We substitute for v^ from (3.10a) into (3.10b) and use the solutions for 


Uq and Uj given above to obtain 


3 2 v 


01 


35' 


8l 5v 01 


V 2 v 
0101 


2miv 0l 2iya + 2„ 2 „ 0l , 

- ®2 fv 01 + i(a +iy 2 ) + y- = 1 V> ? 01 


(3.11) 


Hence, in order to obtain a solution for vqi we require solutions for w m g 
and Vq^. In order to determine w m g we equate coefficients of 6E in 
(2.4d). The resulting equation is 


On substitution of 


U 1 

ia - W m0 



♦ iK l“0 

”„0 + V 01 

u o 

’Si 

3 2 V ( 


3 2 v, 


into (3.12) we obtain 


35 ' 


35 


01 

2 " V 01 


a 2 v, 


01 


a 2 w 


mO 


35 


35 


(3.12) 


and 


01 


from (3.11) and (3.2) respectively 


35 


3 w 

mO ift£ 


iKj/a+2y 2 


w 


mO 


2 M 2 

- 7 w m0 V 0l 


■'V 


(3.13) 
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2ift£ 


211^/ a+2y 2 


+ 3(a+2y 2 ) V 01 V 01 + 3/ -^- v 01 V 01* 


Thus, once we have a solution for Vqj from (3.2), it remains to solve the 
coupled equations (3*11) and (3.13) with the boundary conditions 


v 01* W m0 + 0 as ? * *"• 


(3.14) 


In order to eliminate f(x) from (3.2) we introduce the variable £ 


and let 


g 2 f 

5 =5 i"IT’ 


1 


(3.15) 


so that (3.2) becomes 


9 2 V, 


35 


F- + Si« v oi 

1 


V 01 

~T~ 


(3.16) 


If we look for a solution with Kj of the form 


Qg 2 f>/G 0 x 


K 1 ' K 1 + I ,3/2 * 


g 1 (a+2y 2 )‘ 


(3.17) 


then, using the transformation (3.15), equations (3.11) and (3.13) become, 
respectively, 


a 2 
9 v, 


35 


01 X c 

2 + g l 5 l v 01 
1 


V 01 v 01 2i ^l v 01 2iK 1 /a+2y 2 v Q1 

6 + 3^27? + 


T w m0 V 01 ’ (3,18) 


and 


a 2 

3 w 


35 


mO 

T~ 

1 


iflS l w mO 1C l' ,a+ ^2 ”«0 


a+2y. 


V 


2 V 2 

3 W m0 V 01 


(3.19) 
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+ 


21 ^l''01 V 01 , 
3(a+2y 2 ‘) 


21K 1 /a+2y 2 v 0 jV 01 


VG qX 


In order to simplify (3.16) further we make the transformations 


and 


, .1/3- 

n * (-g^ 


v oi ■ 


(3.20) 

(3.21) 


with the result that (3.16) becomes 



(3.22) 


Note that gj(x), defined by (3.3a) is less than zero for all values of x. 
As noted by HL, (3.22) is a particular form of the second Painleve transcend- 
ent and has been shown by Hastings and Mcleod (1978) to have a solution such 
that 

v n ~ (-n) 1 ^ 2 , n 

and (3.23a,b) 

V ~ /2 Ai(n), n > +». 
n 


A solution for V Q was obtained numerically by using (3.23) and starting in- 
tegrating at n = +“ and integrating to the left with a fourth-order Runge- 
Kutta scheme. These results were used to determine solutions for vg^ and 

w m0* 

For a fixed, real value of 0 we can solve (3.18) and (3.19) and 
determine the complex function K^(x) as x moves downstream. However, we 
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are interested in neutrally stable solutions so we seek real values of ft 
and K^(x). We can eliminate the x-dependence from the coefficients in 
(3 . 18) and (3*19) by making the transformations 


and 


v 


01 



VV <-8l> 2/3 

/a+2y 2 


ft = £2(-g 1 )(a+2y 2 ). 


(3.24) 


(3.25) 


(3.26) 


* a 

Now we can find the constants and £2 so that the flow is neutrally 
stable at the location x where and £2 satisfy (3.25) and (3.26). 
Hence, with the transformations (3.2a), (3.24), (3.25) and (3.26), equations 
(3.18) and (3.19) become, respectively, 


d 2 V 


dn' 


n 710 2iK 1 

£ - (i + — -)TiV p - v p 


V 2 V - w V 
n p 3 w m0 n* 


,2 

d w 


(3. 27a,b) 


dn 


m0 * A 

T ~ 1£iT1W m0 - iK l W mO 


i£2nv V + iK V V . 
mOn 3 pn 3 lpn 


We seek solutions of (3.27), with and £2 real, satisfying 


V W m0 * 0 as 11 * ± “‘ 


As n + +“ (3.27) can be written as 


d 2 V 


2iK, 


dn 


£ _ (1 + ^nVp - -gj- v p = 0, 


(3.28) 
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and 


,2 

d w 


dn 


mO 

2 ~ 


innW tnO ‘ 1K l W mO = 0< 


(3.29a,b) 


Hence, as n +“ we can find two independent solutions for V p and w m Q, 

in terms of the Airy function Ai , which satisfy (3.28). When n + -°> the 

equations for V p and w m Q are 

d 2 V o.n 2 iK i V ,/T 

P _ 2i« v _ 1 , P = _ 2/6 / — 

2 3 n p 3 3 mO ^ ’ 


dn 


,2 

d w 


(3.30a,b) 


dn 


¥- - < lS + 4) " - “iv ■ l5 " / - n v p + • 


The appropriate expansions of (3.30) now take the form 




3/2 


w = h 1 1 / 2 e _<#> I [w + ...] 

w rnO 1 L ™nn J * 


mOO 


where <f> satisfies 

A 

4 ? Si * A ?iQ 

<p + <t> [4 + n] + (8 + in) — 3 - = 0 , 

and we take the two roots of this equation with positive real part to generate 
two independent solutions of (3.30) with Vp, w ffi Q 0, n 

These asymptotic solutions for Vp and w m Q at ^ = ±» were used as 
initial values in the numerical integration scheme used to solve (3.27). 
Equations (3.27) were written as a system of four first order differential 
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equations. This system was solved using a standard fourth order Runge-Kutta 
integration scheme. The integration procedure was started at rj = -°» and 
n = + 00 and continued to n = 0, finding two independent solutions for 
V p and w m Q from each direction. At n = 0 the continuity of a linear 
combination of the independent solutions from each direction produces an 

a a 

eigenvalue problem for and ft. We used a Newton-Rahpson iteration 

A A 

scheme for two variables to find real values for and ft. Using the 

/A A 

above scheme solutions for K^, ft, and w m Q were obtained but we post- 

pone the discussion of these results until the next section. 

A A 

Having found ft and the dimensionless frequency and wavenumber, 

ft and K respectively, can only be found once the HL calculation has been 

performed for a particular curvature distribution x(x). However HL gave 

asymptotic solutions of the free boundary problem for x close to the linear 

neutral position, x = x , and for x a long way downstream of that position 

for curvature distributions which increase as quickly as for x » 1. 

a 

Firstly, we recall that when x x + the shear layers coalesce. Thus 

if we denote by ft T and ft^ the frequencies of the wavy vortex modes 

neutral in the upper and lower layers at x it follows that 


Next suppose that 


* 

ftT* x * x + * 
T 


x M , M > -j for large 


are all given in HL so tht 


(3.31) 

x, the asymptotic forms 
S1T» &1L the values of 


M 

7 


M-l r 
x G. 


’1L 


M 


3M-2 


for a, y^, y 2 and b 
gj are given by 


IT 


(3.32) 
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(We note here in passing 
for x » 1 we obtain 

or 

It follows that as x increases the frequency of the upper layer mode which 
is neutral at x increases whilst that of the lower layer tends to a constant 
value. Thus we can distinguish between the modes as being of high and low 
frequency respectively; this result is entirely consistent with the experi- 
mental results discussed in the next section. 


that gi in the lower layer is positive.) Thus 


0 : M 2M- 12 ft 

~ 9 X G 0’ fl L ~ 3 


MG 0 2M-1 
\ = — X ' 


4. RESULTS AMD DISCUSSION 

The numerical scheme outlined above was used to search for eigenvalues 

A A A A 

(Kpft) in the region > 0, ft > 0. The only eigenvalues located were 

(Kj,a) = (0.97,0.58), (K^ft) = (5.132,1.805). (4.1a,b) 

A A 

It is possible that other eigenvalues exist at higher values of Kp ft 
since for | K 1 I > |ft | > 0 a detailed eigenvalue search was not carried 
out. For any particular incoming boundary layer profile the frequency ft 
and wavenumber K * + e 4 * * 7 are then calculated from (3.5), (3.17), 

(3.25), (3.26) and (4.1). This can only be done once the HL calculation has 
been carried out and the frequency and wavenumber obtained in this way will of 
course be dependent on x. The resulting expressions for ft and K should 
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be interpreted as the frequency and wavenumber which are neutrally stable at 
x. Alternatively we could invert the equation 

ft = ft(x,e ) , 

to find the downstream location where the wavy vortex is neutrally stable. 
If ft is held fixed at the neutral value at x = x then the wavenumber 
K becomes complex for x * x so the wavy vortex mode undergoes spatial 
amplification or decay away from the neutral location. 

We do not repeat the HI, calculation in order to obtain specific values 
for K, ft for a particular boundary layer flow. We believe that the major 
result of this paper is that the large amplitude states of HL are unstable in 
the thin shear layers which trap the vortices. The only available experi- 
mental results which give detailed results on the wavy mode structure do not 
give sufficient detail about the unperturbed boundary layer to enable us to 
calculate the relevant values of K and ft, so below we discuss only the 
qualitative agreement between our theory and these experimental results. How- 
ever, before we compare our results with experiments we shall first describe 
the eigenfunctions appropriate to (4.1). 

These functions are shown in Figures 2, 3 and we point out their oscilla- 
tory nature for negative values of q. We note that q + -» corresponds 
to moving from IIa,b into the core region I. It is of course possible that 
other "lower order" modes with a less oscillatory nature might exist but, as 
stated already, only those corresponding to (4.1) were found. 

Now let us turn to the physical implications and experimental relevance 
of our calculation. We first stress that the instability mechanism which we 
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have described in detail for just the upper shear layer can also occur in the 
lower layer* The modes of instability of the shear layers are independent be- 
cause they decay exponentially away from the center of the layers* Thus if we 
consider the frequency of the imposed wavy mode to be fixed then the layers 
will breakdown in the manner described at different downstream locations. 
Since the downstream velocity component of the basic state is largest in the 
upper shear layer it is to be expected that this layer will be the first to 
become unstable. Furthermore, since the wavenumbers K(x) appropriate to a 
fixed frequency disturbances will be different, the modes propagates down- 
stream with different wavespeeds. In fact initially, by which we mean close 
to the linear neutral position, the layers IIa,b coalesce so that if breakdown 
occurs close to this point then the structure in these layers will be very 
similar. Further downstream the upper layer moves into the free stream and 
the downstream velocity component tends to the free stream speed. The lower 
layer however approaches the boundary so that the fluid velocity there tends 
to zero. Thus, it follows that if the stationary Gortler vortices develop 
over a sufficiently long interval before breakdown then the upper layer wavy 
mode moves downstream with the free stream speed whilst the lower one has a 
much smaller propagation speed. It is interesting to note that in the 
apparently closely related Taylor vortex problem the corresponding breakdown 
is due to a single wavy mode whose presence is felt throughout the flow. 

There have been many experimental investigations of the secondary insta- 
«« 

bility of Gortler vortices; the reader is referred to the papers by for 
example Bippes and Gortler (1972), Wortmann (1969), and Bippes (1978) who all 
described the secondary mode as being locally periodic in x and t. However 
more recently Kohama (1987) and Peerhossaini and Wesfreid (1987) have given 



-23- 


more details of the flow structure which exists when breakdown occurs. We 
first discuss the results of Peerhossaini and Wesfreid. 

The boundary layer investigated by the latter authors was in the concave 
section of a curved channel. They found that when the secondary mode first 
appeared it was confined to a region at the top of the vortices. We interpret 
this as being due to the instability mechanism we have described being first 
operational in region Ila. Further downstream they reported a similar insta- 
bility but this was localized near the wall. We interpret this instability as 
being due to the wavy vortex instability of region lib. Both instabilities 
observed by Peerhossaini and Wesfreid had wavy vortex boundaries in the down- 
stream direction, this is entirely consistent with the breakdown mechanism we 
have described in Section 3. 

The experiments of Kohama were performed on the NASA laminar flow wing 
discussed by for example Pfenninger, Reed, and Dagenhart (1980). Kohama 
described only the breakdown in the upper part of the boundary layer but gave 
measurements of the wavespeeds at different downstream locations. In the lam- 
inar flow region of the wing the wavespeed was found to be about 0.45 of the 
free stream speed but this factor became about 0.99 when the flow was fully 
turbulent. Such a variation of wavespeed is predicted qualitatively by our 
theory since the layer Ila migrates from being somewhere in the middle of the 

boundary layer (more precisely at the position where Rayleigh's criterion for 
• • 

the pre-Gortler flow is most violated) to the edge as the vortices become ful- 
ly nonlinear. The migration of this layer we believe accounts for the varia- 
tion of wavespeed given by Kohama. 

Finally we note that there are other instability mechanisms which could 
account for the appearance of three-dimensionality and time-dependence in the 
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later stages of Gortler vortex development. The most likely other types of 
disturbances would be Rayleigh instabilities associated with the spanwise 
locally inflexional velocity profiles which certainly develop as the vortices 
develop in a steady manner. Secondly, there exists the possibility that 
Tollmien-Schlichting waves might cause the flow to become three-dimensional. 
However, it is not clear that these modes would lead to the wavy vortex 
boundaries which seem to be always observed when breakdown occurs. Never- 
theless, it is likely that in some situations the Rayleigh and Tollmein- 
Schlichting modes might be important in the later stages of the transition 
process . 
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